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ABSTRACT 


It  is  the  purpose  of  this  paper  to  review  the  recent  history  of  the  development  of  py¬ 
rolytic  graphite  and  its  thermal,  aerothermochemical,  and  structural  properties 
with  regard  to  its  application  to  hypersonic  vehicles. 

Pyrolytic  graphite  is  not  really  new  since  its  discovery  dates  back  to  Thomas  Edison, 
however,  this  class  of  refractory  materials  has  been  "rediscovered"  recently  because 
it  possesses  certain  desirable  properties  at  extremely  high  operating  temperatures, 
including  a  high  strength  to  weight  ratio,  low  oxidation  rate,  and  a  low  transverse 
coefficient  of  thermal  conductivity. 

Consequently,  a  thorough  investigation  of  the  properties  of  pyrolytic  graphite  is  being 
carried  out.  In  addition,  theoretical  models  are  being  developed  so  that  its  aero¬ 
thermochemical  and  aerothermostructural  behavior  can  be  predicted  with  confidence 
for  arbitrary  hypersonic  missions.  In  this  paper,  the  currently  available  property 
data,  and  analytical  and  numerical  techniques  are  presented  and  discussed  critically, 
and  areas  for  future  research  and  development  are  delineated. 
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PYROLYTIC  GRAPHITE  - 
A  STATUS  REPORT 

A.  M.  Garber,  E.  J.  Nolan  and  S.  M.  Scala 
General  Electric  Company 
Missile  and  Space  Division 
Philadelphia  1,  Pennsylvania 


INTRODUCTION 

One  of  the  newer  refractory  materials  which  has  received  a  certain  amount  of 
emphasis  for  hypersonic  applications  during  the  past  few  years  is  pyrolytic  graphite. 

Surprisingly,  pyrolytic  graphite  is  not  really  new  since  its  discovery  dates  back 
to  Thomas  Edison  (1).  One  of  the  interesting  features  of  pyrolytic  graphite  which 
makes  it  attractive  for  hypersonic  vehicles  is  the  high  degree  of  anisotropy  of  this 
material  in  both  its  physical  and  mechanical  properties  and  its  low  rate  of  oxidation 
(ablation). 

Indeed,  several  years  ago,  its  use  was  considered  in  one  of  the  active  missile 
programs,  but  because  of  certain  logistics  problems  that  arose,  it  was  not  chosen 
for  flight  hardware.  Intensive  materials  efforts  during  this  particular  program  gave 
rise  to  the  first  full  scale  manufacturing  effort  based  on  the  vapor  deposition  of  py¬ 
rolytic  graphite.  Prior  to  this  time,  only  small  pieces  had  been  made  on  a  labora¬ 
tory  scale,  by  Edison  in  his  incandescent  lamp  work,  and  more  recently  by  several 
oilier  investigators. 

Pyrolytic  graphite  may  be  produced  by  thermally  decomposing  a  hydrocarbon 
gas,  such  as  methane,  and  depositing  the  carbon  on  a  hot  surface,  under  carefully 
controlled  conditions.  One  method  now  in  use  consists  of  heating  the  hydrocarbon 
gas  to  temperatures  ranging  up  to  5000°R,  in  a  vacuum  furnace,  see  Fig.  (1).  Dur¬ 
ing  the  decomposition  of  the  gas,  gaseous  carbon  deposits  on  the  surface  of  a  pre¬ 
pared  mandrel  which  is  usually  commercial  graphite.  The  carbon  atoms  formed  dur¬ 
ing  the  cracking  process,  then  condense  on  the  substrate  in  an  orderly  fashion,  layer 
upon  layer,  producing  an  ordered  structure  which  has  a  higher  strength  to  weight 
ratio  than  commercial  graphite,  see  Fig.  (2),  and  hence  is  potentially  more  attractive 
for  aerospace  applications.  Some  typical  specimens  of  pyrolytic  graphite  are  shown 
in  Fig.  (3). 

The  separate  hexagonal  layer  planes  of  pyrolytic  graphite  are  all  parallel  to  the 
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Fig.  2.  Strength  to  Weight  Ratio  vs.  Temperature 
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basal  (deposition)  plane  (a-b  plane),  but  do  not  have  the  regular  appearance  of  a  honey¬ 
comb  with  respect  to  each  other.  Microscopic  examination  shows  that  pyrolytic  gra¬ 
phite  is  a  highly  oriented  crystalline  structure  resembling  densely  packed  stalagmites, 
see  Fig.  (4). 

In  the  case  of  commercial  graphite,  the  layer  planes  have  a  random  orientation 
for  the  different  crystallites;  there  are  approximately  two  crystals  oriented  normal 
to  the  surface  ("c"  axis),  for  every  crystal  oriented  parallel  to  the  surface  (a-b  plane). 
In  contrast,  the  crystals  of  pyrolytic  graphite  are  highly  oriented,  so  that  the  number 
oriented  in  the  "c"  direction  is  two  or  three  orders  of  magnitude  for  every  one  in  the 
a-b  orientation.  The  high  degree  of  anisotropy  in  the  structure  of  pyrolytic  graphite 
results  in  a  high  degree  of  anisotropy  in  all  of  its  thermal  properties,  (except  the  spe¬ 
cific  heat). 

During  the  early  manufacturing  phases  of  the  full  scale  hardware  of  pyrolytic 
graphite  in  the  missile  program  mentioned  above,  certain  problems  arose  that  were 
not  evident  during  the  laboratory  effort,  resulting  in  what  appeared  to  be  random 
cracldng,  delamination  and  spallation  of  the  freestanding  piece  as  manufactured.  It 
was  observed  further  that  hardware  which  appeared  to  be  of  sound  quality  would  often 
fail  upon  exposure  to  heating  in  arcs  and  rocket  engine  exhausts.  The  nature  of  these 
failures  was  similar  to  those  occurring  in  the  freestanding  manufactured  piece.  An 
exhaustive  investigation  into  the  causes  of  these  failures  was  undertaken  by  the  manu¬ 
facturers  involved  in  this  missile  program  and  within  six  months,  a  number  of  inter¬ 
esting  results  were  obtained. 

It  was  found  for  example,  that  nodule  formation,  which  was  a  source  of  stress 
concentration,  and  therefore  potential  failure,  could  be  minimized  by  improving  cer¬ 
tain  process  parameters  and  the  housekeeping  procedures  around  the  process  equip¬ 
ment.  The  question  as  to  why  delamination  and  cracking  of  the  manufactured  pyrolytic 
graphite  can  occur,  involves  a  thorough  understanding  of  the  behavior  of  brittle  aniso¬ 
tropic  materials,  the  technique  of  deposition  and  the  subsequent  stresses  which  arise 
during  the  deposition  process  and  the  cooling  down  period.  This  stress  interplay  and 
its  importance  was  recognized,  (2)  during  this  investigation  period.  Techniques  for 
computing  the  magnitude  of  these  stresses  in  an  anisotropic  material  could  not  be 
developed  in  such  a  short  period  of  time.  However,  enough  of  the  overall  problem, 
including  proper  design  of  the  specimen  under  test,  was  sufficiently  well  understood  to 
permit  some  thirteen  straight  successful  experiments  to  be  performed  in  a  15"  dia¬ 
meter  Vanguard  rocket  engine  exhaust.  It  is  noteworthy  that  full  scale  experimental 
hardware  was  used  during  these  tests. 

During  and  subsequent  to  the  aforementioned  pyrolytic  graphite  development 
effort,  the  investigation  of  analytical  techniques  for  the  prediction  of  mass  transfer, 
heat  transfer  and  the  structural  response  of  pyrolytic  graphite  during  atmospheric 
flight  at  hypersonic  speeds  was  initiated  at  the  General  Electric  Company.  Several 
reports  of  different  phases  of  this  effort  have  already  appeared  in  the  literature,  see 
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(3),  (4),  (5)  and  (6).  It  has  been  found,  for  example,  that  by  correct  mandrel  design 
and  control  of  certain  variables  during  deposition,  freestanding  pieces  can  be  prepared 
that  would  experience  thermal  stresses  during  flight  which  are  lower  in  magnitude  than 
those  obtained  during  manufacture  and  the  subsequent  cool-down  process.  It  has  been 
found  also  that  a  change  in  the  degree  of  anisotropy  can  likewise  change  the  thermal 
stress  distribution  produced  during  flight,  depending  on  the  level  and  direction  of  the 
anisotropy. 

These  results  will  be  reviewed  in  more  detail  in  a  later  section. 

CHEMICAL  KINETICS 

In  considering  the  chemical  aspects  of  the  performance  of  pyrolytic  graphite,  one 
quickly  notes  a  certain  paradoxical  situation.  That  is,  although  the  chemical  reactions 
between  carbon  and  oxygen  have  been  studied  for  well  over  one  hundred  years,  never¬ 
theless  an  incomplete  understanding  of  the  mechanism  of  mass  loss  exists  even  for 
ordinary  graphite.  Although  it  is  known  for  a  certainty  that  in  the  low  surface  tempera¬ 
ture  regime  where  the  oxidation  process  is  rate  controlled,  the  volatile  reaction  pro¬ 
ducts  CO2  and  CO  are  produced,  it  has  not  yet  been  established  as  to  what  the  actual 
sequence  of  steps  is,  nor  has  the  slowest,  and  hence  rate  determining,  step  been  deter¬ 
mined.  Consequently,  it  is  both  useful  and  necessary  to  introduce  the  following  empiri¬ 
cal  Arrhenius  representation  for  the  oxidation  rate  in  the  reaction  rate  controlled  oxi¬ 
dation  regime 


(1) 


where  k  is  the  specific  reaction  rate 
oxygen  near  the  surface,  and  n  is  the  order  of  the 


"  W.* 


the  partial  pressure  of  molecular 


reaction. 


An  examination  of  the  available  experimental  data  0,  (8),  indicates  that  although 
n  is  bounded,  there  is  some  disagreement  as  to  the  precise  magnitude  of  n.  Further, 
upon  writing  the  specific  reactivity  in  the  form  (9), 


,  .  -E/RTW 

k  =  k  e  w  (2’ 

o  '  ‘ 

it  is  seen  that  in  order  to  specify  the  chemical  kinetic  behavior  of  pyrolytic  graphite, 
one  requires  a  knowledge  of  E,  the  activation  energy  and  k0,  the  effective  collision 
frequency. 


8 


Upon  making  use  of  the  data  appearing  in  Ref.  (9),  one  finds  that  for  graphitic 
materials,  the  following  bounds  on  the  chemical  kinetic  properties  apply:  (see  Fig.  (5)) 


0  <  n  <  1.0 

(3) 

8  <  E  <  60  Kcal/ mole 

(4) 

1<  kQ<  109  lb. /ft. 2  sec. 

(5) 

This  extremely  wide  range  in  chemical  kinetic  data  for  graphite  may  be  attributed  to 
the  fact  that  these  properties  depend  critically  on  how  the  material  was  manufactured. 
Thus,  since  there  are  at  least  as  many  different  grades  of  graphite  as  there  are  manu¬ 
facturers  (with  different  quality  control  procedures),  it  is  unlikely  that  any  two  investi¬ 
gators  have  tested  the  same  material. 

At  the  moment,  the  only  useful  data  on  the  oxidation  rate  of  pyrolytic  graphite 
is  due  to  Horton  (10).  However,  as  has  been  pointed  out  by  Nolan  and  Scala  (4),  there 
is  sufficient  scatter  in  his  data  to  question  his  optimistic  (extremely  low  mass  loss) 
interpretation  of  the  data.  Consequently,  at  this  moment  we  can  quote  two  sets  (4), 
(optimistic  and  conservative)  of  chemical  kinetic  data,  based  on  one  set  (10),  of  experi¬ 
mental  measurements  for  the  oxidation  of  pyrolytic  graphite. 


E  =  32. 8  Kcal/mole 
k  =  3.  07  lb.  /ft.  2  sec. 


(6a) 


E  =  47.  5  Kcal/mole 

k  =  1.28  x  10^  lb.  /ft. 2  sec. 
o 


(6b) 


It  is  noted  that  the  value  of  n  to  be  used  with  the  data  of  Eq.  (6a)  or  (6b)  is  n  =  1/2. 
For  the  case  of  n  =  1/2,  and  for  low  surface  temperatures,  in  the  rate  controlled  oxi¬ 
dation  regime,  since  the  pressure  gradient  at  the  surface  is  negligible,  one  can  write 
for  the  partial  pressure  of  oxygen  at  the  surface, 

(P0  “  (\)  Pe 

w  w 

and  hence,  there  follows  immediately: 
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Rate  Controlled  Oxidation  Data,  (P  =  0.  1  atm. 


where,  in  first  approximation, 


w 


It  is  expected  that  Eq.  (8)  will  apply  even  when  the  free  stream  air  is  dissociated  pro¬ 
vided  that  the  gas  phase  reaction  rates  are  sufficiently  fast  so  that  the  atomic  oxygen 
has  recombined  to  the  molecular  state  before  colliding  with  the  surface.  If  the  gas 
phase  reaction  rates  are  either  "frozen",  or  so  slow  that  oxygen  atom  recombination 
is  not  complete  in  the  gas  phase,  then  the  rate  constants  given  in  Eq.  (6)  are  not  appli¬ 
cable,  since  atomic  and  molecular  oxygen  react  with  carbon  at  different  rates  (11). 

THERMAL  AND  STRUCTURAL  PROPERTIES 

Even  though  pyrolytic  graphite  is  not  a  new  material,  a  significant  attempt  to 
measure  its  properties  has  only  been  made  recently,  (12)  through  (22),  and  a  review 
of  information  available  since  1959  on  certain  properties  of  pyrolytic  graphite  has  re¬ 
cently  been  made,  (23).  For  data  to  be  reliable,  and  to  insure  proper  deductions  from 
theoretical  calculations,  it  is  essential  that  the  experimental  data  be  obtained  for  pyrolytic 
graphite  that  has  been  manufactured  under  controlled  conditions.  In  particular,  all  of 
the  properties  should  be  determined  for  specimens  from  the  same  manufactured  batch. 

However,  in  practice,  it  appears  that  the  acquisition  of  properties  of  a  particular 
identifiable  pyrolytic  graphitic  material  is  difficult.  Because  of  manufacturing  and  phy¬ 
sical  variables  that  occur  during  deposition,  such  as  grain  size,  pressure  changes  and 
changes  in  flow  rate,  one  may  expect  that  one  type  of  pyrolytic  graphite  manufactured 
in  one  facility  need  not  necessarily  be  the  same  as  that  from  another.  As  a  consequence 
of  these  process  variations,  substantial  variations  in  the  property  values  have  been  re¬ 
ported  in  the  literature.  Unfortunately,  in  some  instances,  the  data  reported  in  the 
literature  does  not  identify  the  material  from  which  the  data  was  obtained. 

In  reviewing  the  literature,  it  has  been  found  also  that  the  details  of  the  experi¬ 
mental  procedure  have  not  been  given  adequately.  For  example,  the  preparation  of 
samples  may  introduce  delaminations,  which  can  affect  the  strength. 

It  should  be  noted  that  it  is  difficult  to  measure  stresses  and  strains  in  pyrolytic 
graphite  because  of  its  brittle  nature  at  relatively  low  temperatures.  In  making  these 
measurements  small  bending  stresses  can  appreciably  alter  measured  fracture  stresses. 
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When  making  measurements  at  elevated  temperatures  difficulties  are  introduced 
because  of  anisotropic  expansion,  and  further  difficulties  arise  because  of  the  lack 
of  a  good  technique  in  accurately  measuring  high  surface  temperatures  in  a  range  from 
4500°R  to  5500°R.  It  should  also  be  mentioned  that  certain  qualifications  are  necessary 
prior  to  the  application  of  data  obtained  on  flat-wall  machined  specimens  to  curved 
shapes  such  as  might  be  found  on  hypersonic  vehicles. 

It  has  been  observed  (23)  that  machining  can  introduce  surface  defects  which  might 
induce  fracture  failures  far  sooner  than  would  be  the  case  with  smooth  undamaged  sur¬ 
faces.  Of  course,  this  does  not  apply  to  such  values  as  the  elastic  modulus,  thermal 
conductivity  and  coefficient  of  expansion.  When  measuring  the  high  temperature  pro¬ 
perties  (above  3C00  R),  one  should  also  take  into  account  the  probability  that  pyrolytic 
graphite  becomes  translucent  and  that  phonon  transport  plays  a  significant  role  in  de¬ 
fining  the  apparent  thermal  conductivity. 

By  far,  the  largest  amount  of  data  available  has  been  published  as  a  Lockheed 
Missiles  and  Space  Company  report  (20).  Unfortunately,  however,  the  data  presented 
is  not  readily  identifiable  with  a  well-characterized  material,  and  hence  its  usefulness 
is  severely  limited.  Insofar  as  data  reported  in  the  literature  are  concerned,  the  re¬ 
ports  issued  by  the  Raytheon  Company  (21)  and  Aerojet  General  (22)  are  quite  complete. 
Their  evaluation  of  the  experimental  difficulties  involved  ir  making  pyrolytic  graphite 
measurements  and  the  necessity  of  correlating  measured  data  to  a  well-defined  material 
show  an  appreciation  for  the  experimental  techniques  and  manufacturing  difficulties  asso¬ 
ciated  with  anisotropic  materials.  However,  this  is  not  to  imply  that  the  materials  used 
in  their  experiments  were  always  satisfactory.  Indeed,  the  reports  (22)  state  that  spe¬ 
cimens  were  received,  and  presumably  used,  in  which  defects  such  as  chips,  partial 
delaminations  and  large  nodules  were  visible. 

The  only  data  reported  which  considered  the  effect  of  structure  on  properties  were 
found  in  a  summary  report  issued  by  Raytheon,  (21).  A  microsti’uctural  classification 
technique  was  used  to  identify  the  material.  Both  mechanical  and  thermal  properties 
were  obtained  for  several  different  types  of  microstructures  which  varied  from  surface 
nucleated  to  highly  regenerative  material.  The  majority  of  the  data  is  based  on  high¬ 
ly  regenerative  material,  but  enough  work  was  performed  on  a  surface  nucleated 
structure  to  allow  comparison.  These  two  structures  differ  considerably  from  each 
other  and  hence  may  be  taken  as  approximate  bounds  when  considering  the  property 
variations  of  pyrolytic  graphite. 

The  coefficient  of  thermal  expansion  for  surface  nucleated  and  the  highly  regenera¬ 
tive  continuously  nucleated  materials  in  both  the  "a"  and  "c"  directions  is  shown  in 
Figs.  (6)  and  (7).  The  regenerative  material  which  has  more  contributions  from  the  "c" 
direction,  shows  a  significantly  higher  coefficient  of  expansion.  On  the  other  hand,  the 
"c"  direction  coefficient  of  expansion  (Fig.  7)  seems  to  be  independent  of  microstructure. 


The  thermal  conductivity  coefficient  for  both  the  "a"  and  "c"  directions  is  shown  in 
Figs.  (8)  and  (9).  The  regenerative  material  has  a  slightly  higher  conductivity,  which 
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Fig.  7.  Coefficient  of  Thermal  Expansion  of  Pyrolytic  Graphite  in  "C"  Direction 

for  Various  Microstructures 
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may  or  may  not  be  extrapolated  to  higher  temperatures  because  of  the  increasing  role 
of  phonon  transport.  The  difference,  however,  is  not  large.  It  should  be  pointed  out 
that  heat  treatment  at  2900°C  for  two  hours  can  increase  the  "a"  direction  conductivity 
by  about  a  factor  of  three,  (21). 

Dynamic  modulus  data  (21)  are  shown  in  Fig.  (9).  The  data  are  interesting  from 
several  respects.  An  anomalous  increase  in  modulus  in  a  range  of  temperatures  from 
250°  to  700°F  is  shown  together  with  the  expected  linear  decrease  in  modulus  with  an 
increase  in  temperature  above  700°F.  In  the  absence  of  relaxation,  this  linear  decrease 
is  to  be  expected  for  normal  polycrystalline  materials.  Fig.  (10)  shows  that  the  mo¬ 
dulus  decreases  linearly  from  17500  to  2250OF;  but  between  1000O  and  1750OF  shows 
a  departure  from  linearity.  This  type  of  behavior  has  often  been  noted  for  polycrystal¬ 
line  materials  and  is  attributed  to  grain  boundary  plastic  flow.  In  pyrolytic  graphite, 
the  behavior  may  be  due  to  microscopic  viscous  flow  between  conical  crystallite  bound¬ 
aries.  If  this  is  indeed  the  case,  some  relaxation  of  residual  stresses  may  be  expected 
in  this  temperature  range. 

The  dynamic  modulus  of  the  surface  nucleated  material  is  slightly  higher  than  the 
regenerative  material  throughout  the  temperature  range  75°  to  2250°F.  The  static 
room  temperature  moduli,  also  shown  on  Fig.  (10)  are  somewhat  higher  but  bear  the 
same  relation  to  one  another,  as  do  the  dynamic  moduli. 

The  difference  in  magnitude  between  static  and  dynamic  values  was  not  discussed 
in  Ref.  (21),  but  may  well  be  due  to  using  the  resonant  frequency  obtained  experiment¬ 
ally  and  computing  the  modulus  on  the  basis  of  isotropic  elastic  theory.  In  general, 
static  and  dynamic  moduli  for  isotropic  polycrystalline  metals,  in  the  absence  of  relaxa¬ 
tion,  are  equal  (24). 


It  is  not  realistic  to  use  the  temperature  dependence  of  the  dynamic  modulus  and 
apply  it  to  the  static  values  because  of  the  (presumed)  conical  boundary  flow  between 
750°  and  1750°F.  This  behavior,  whatever  the  cause,  leads  to  a  45%  drop  in  the  dy¬ 
namic  modulus  between  75°F  and  2250°F.  On  the  other  hand,  recent  data  obtained  by 
GE  show  that  the  static  modulus  of  a  highly  regenerative  material  is  essentially  con¬ 
stant  from  75°F  to  3000°F  changing  only  from  3.  5  to  3.  4  x  10®psi.  The  cause  of  the 
temperature  variation  of  the  dynamic  modulus,  and  associated  variation  of  the  internal 
friction,  is  not  fully  understood  but  should  be  studied  further. 


The  values  reported  in  Ref.  (21)  for  the  static  elastic  modulus  of  the  regenerative 
structure  are: 


Static  Modulus  of 
Regenerative  Material 


Tension,  psi. 


Compression,  psi. 


"a"  Direction 

4.0-4.4x106 

4.  lx  106 

"c"  Direction 

1.  6x  106 

1.6x106 

These  values  are  slightly  higher  than  the  foregoing  data,  but  indicate  that  the  tensile 
and  compressive  moduli  are  similar.  The  difference  between  the  data  may  well  be 
due  to  different  manufacturing  techniques. 

Kotlensky  and  co-workers  at  JPL  have  published  several  reports  dealing  with  PG 
over  the  past  several  years  (25)  through  (28).  The  data  include  only  tensile  tests,  and 
are  shown  in  Table  (1).  The  data  are  plotted  in  Fig.  (11),  together  with  50%  confidence 
limits  (20).  The  Jet  Propulsion  Laboratory  data  seem  to  conform  to  the  50%  curve. 
These  data  also  show  mechanical  properties  of  pyrolytic  boron  nitride.  It  will  be  noted 
that  the  low  temperature  strength  appears  to  be  lower  than  that  of  pyrolytic  graphite. 

Work  has  also  been  performed  recently  (29)  at  General  Electric  -  MSD.  Data  ob¬ 
tained  (but  not  presented  here),  include  both  continuously  nucleated  and  surface  nucle¬ 
ated  material.  A  report  is  currently  in  preparation  and  stress-strain  curves  will  be 
available.  Other  data  obtained  as  part  of  the  aforementioned  work  are  the  coefficient 
of  thermal  expansion,  shear  strength  in  the  a-b  plane,  and  three  and  four  point  loaded 
bend  tests,  all  measured  at  temperatures  up  to  5000OF.  The  new  data  for  pyrolytic 
graphite  shows  appreciably  higher  strength  than  both  the  Lockheed  mid  JPL  data.  As 
a  consequence  of  this  study,  these  investigators  are  lead  to  conclude  that  "pyrolytic 
graphite  as  it  is  currently  made  and  tested,  is  a  considerably  stronger  and  more  re¬ 
liable  material  than  was  the  case  earlier  in  its  development".  This  may  be  attributed 
to  both  improved  process  control  and  a  better  understanding  of  the  effects  of  surface 
damage  on  the  properties  of  pyrolytic  graphite.  Care  in  selection  of  representative 
material,  inspection  of  specimens  before  testing  mid  of  fracture  surfaces  after  testing 
must  become  part  of  the  routine  test  procedure  to  insure  that  quoted  results  corres¬ 
pond  to  the  mode  of  failure  actually  produced  mid  that  they  do  not  reflect  the  effects  of 
notches,  nodules,  delaminations  or  machine  damage. 

It  is  further  believed  by  these  investigators  that  the  maximum  inherent  strength  of 
pyrolytic  graphite  has  yet  to  be  achieved  in  test  specimens.  Studies  of  the  basic  pro¬ 
cesses  by  which  fracture  is  initiated  and  propagated  have  lead  to  the  conclusion  that 
cutting  and  machining  of  test  specimens  leads  inevitably  to  the  presence  of  stress  con¬ 
centration  on  all  surfaces.  "The  strength  as  measured  is  thus  not  indicative  of  the  in¬ 
herent  strength  of  pyrolytic  graphite,  but  rather  of  a  specimen  which  has  already  been 
damaged".  Improved  machining  and  finishing  techniques  are  currently  under  study  as 
well  as  methods  of  restoring  surfaces  to  a  relatively  undamaged  status.  However,  there 
is  a  need  for  developing  a  general  design  and  fabrication  technology  for  pyrolytic  gra¬ 
phite  as  well  as  for  the  newer  vapor-deposited  refractory  miisotropic  materials  which 
are  currently  undergoing  exploratory  evaluation.  The  demand  for  some  of  their  unique 
properties  is  pressing  enough  to  carry  their  study  beyond  the  preparative  stage  to  a 
level  where  they  are  practical  materials  of  construction. 

A  recent  review  article  (23)  has  summarized  the  difference  in  behavior  between 
pyrolytic  graphite  (PG),  boron  pyrolytic  graphite  (BPG),  and  boron  nitride  (BN),  with 
respect  to  thermal  conductivity.  These  data  arc  shown  in  Fig.  (12).  In  addition,  the 
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THERMAL  CONDUCTIVITY,  BTU  -  FT  /  FT2-  °F  -  HR 


Fig.  12.  Thermal  Conductivity  of  Pyrolytic  Materials 
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specific  heat  of  PG  is  shown  in  Fig.  (13).  Because  the  PG,  BPG,  and  BN  are  not  suf¬ 
ficiently  characterized,  these  data  should  be  used  for  comparative  purposes  alone. 

OXIDATION  REGIMES 

If  one  makes  use  of  the  Arrhenius  representation  for  the  rate  controlled  oxidation 
regime  as  discussed  earlier,  it  is  easily  seen  that  as  the  surface  temperature  rises 
above  2000°R,  the  rate  of  mass  loss  increases  exponentially  and  soon  depletes  the  local 
oxygen  supply  which  is  available  in  the  hypersonic  boundary  layer.  Consequently,  the 
slowest  step  is  no  longer  the  specific  chemical  reaction  at  the  surface,  but  rather  the 
rate  at  which  oxygen  bearing  species  are  transported  to  the  surface  by  diffusion  and 
convection. 

Thus,  one  can  distinguish  several  different  oxidation  regimes.  At  the  lowest  sur¬ 
face  temperatures,  one  has  the  reaction  rate  controlled  regime.  But,  as  the  surface 
temperature  rises,  transition  to  a  diffusion  controlled  oxidation  regime  must  occur, 
in  which  the  mass  loss  is  relatively  independent  of  surface  temperature.  Finally,  at 
sufficiently  high  surface  temperatures  (Tw  >  5000°R),  the  rate  of  sublimation  of  chemi¬ 
cal  species  such  as  C(g)  and  C3(g)  exceeds  the  oxidation  rate  and  the  process  is  then 
one  of  sublimation.  For  graphite,  the  term  "sublimation"  should,  of  course,  be  re¬ 
served  for  the  regime  of  extremely  high  surface  temperatures.  In  Fig.  (14),  we  have 
shown  a  representation  of  these  different  regimes.  Note  that,  eventually,  at  suffi¬ 
ciently  high  surface  temperatures,  the  vapor  pressure  can  exceed  the  stagnation  pres¬ 
sure  and  the  boundary  layer  approximation  ^p  _  ^  is  no  longer  valid. 

Sy 

An  exact  solution  for  the  mass  loss  has  recently  been  obtained  by  Scnla  (9),  in  the 
diffusion  controlled  regime.  This  may  be  expressed: 


and  is  shown  in  Fig.  (15),  It  is  noted  that  by  definition,  barring  mass  loss  by  mechani¬ 
cal  spallation,  this  equation  for  mass  loss  applies  equally  well  to  pyrolytic  graphite,  ATJ, 
and  in  fact,  all  grades  of  carbon  and  graphite. 

The  corresponding  heat  transfer  rate  to  the  pyrolytic  graphite  in  the  diffusion  con¬ 
trolled  oxidation  regime  is  given  by  (9): 
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Fig„  14,  Mass  Transfer  Regimes  for  Pyrolytic  Graphite 
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Fig„  15*  Mass  Transfer  in  Rate  Controlled  and  Diffusion  Controlled  Oxidation  Regimes 


Here  Rg  is  the  nose  radius,  Pe  is  the  local  pressure  over  the  object,  A  is  the  angle 
of  yaw  and  6  =  1  for  an  axi-symmetric  body,  while  6  =  0  for  a  two-dimensional  body. 

A  reasonable  representation  for  f ( A )  is  given  by: 

f  (A)  =  cos3//2  A  (13). 

For  the  transition  regime,  Scala  has  shown  that  one  may  obtain  an  approximate  so¬ 
lution  for  the  mass  transfer  in  the  form: 


1/2 


1  +  -  1 


"V2  <V 


(14) 


while  Eq.  (12)  for  the  heat  transfer  rate  applies  to  the  transition  and  rate  controlled 
regimes  as  well. 


Exact  solutions  for  the  sublimation  regime  have  just  recently  been  obtained  by 
Scala  and  Gilbert  (30).  In  this  regime,  one  finds  that  in  order  to  account  properly  for 
the  dominant  aerothermochemical  effects,  a  model  is  required  in  which  one  includes 
the  relative  motion  of  nine  chemical  species  including  0,  02,  N,  N2,  CO,  C02,  CN, 

C(g),  and  C;j(g).  The  first  four  species  are  the  four  major  components  of  dissociated 
iar;  carbon  monoxide,  carbon  dioxide  and  cyanogen  are  combustion  products,  where¬ 
as  C3  and  C  are  the  most  important  gaseous  vaporization  products.  In  Fig.  (16),  the 
concentrations  of  the  predominant  species  present  in  the  gas  adjacent  to  a  graphite 
surface,  are  shown  as  a  function  of  surface  temperature,  for  a  pressure  of  one  atmosphere. 

When  the  appropriate  boundary  layer  equations  are  solved,  one  finds,  as  expected, 
that  the  mass  rate  of  sublimation  increases  exponentially  with  surface  temperature. 

The  effect  of  an  increase  in  stagnation  pressure  is  to  suppress  the  rate  of  mass  loss 
due  to  vaporization  and  the  sublimation  regime  is  shifted  to  higher  surface  tempera¬ 
tures.  Typical  results  (30),  are  shown  in  Fig.  (17)  for  stagnation  pressures  of  10“2  and 
102  atmospheres.  Here  the  total  mass  has  been  normalized  by  mD.  There  is  a  corres¬ 
ponding  significant  reduction  in  the  heat  conducted  into  the  surface,  since  the  mass 
transfer  process  in  the  sublimation  regime  is  highly  endothermic.  Complete  details 
are  given  in  Ref.  (30). 


TRANSIENT  THERMAL  RESPONSE 

The  physicochemical  behavior  of  pyrolytic  graphite  may  be  obtained  by  the  solution 
of  the  energy  equation  applied  to  a  solid  whose  surface  is  undergoing  an  oxidation  re¬ 
action. 
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Graphite  Surface,  (P  =  1  atm. 


1 


(15) 


P'pBr  +v-Vt)  =  -V.Q-  Ew.h. 

where  w.  is  the  chemical  source  term  and  Q  is  the  heat  flux  vector.  Since  the  chemi¬ 
cal  reactions  are  confined  to  the  surface,  the  heat  release  due  to  chemical  reaction  may 
be  consequently  assigned  as  a  boundary  condition,  because  there  are  no  chemical  oxida¬ 
tion  reactions  occurring  within  the  body  of  the  solid.  As  a  consequence,  the  chemical 
source  term  wj  need  not  appear  in  the  fundamental  equation.  If  one  selects  a  coordi¬ 
nate  system  which  remains  stationary  with  respect  to  the  moving  surface  rather  than 
a  coordinate  system  which  moves  with  the  moving  surface,  the  convective  term 
v  •  V  T  is  not  required  and  the  surface  oxidation  appears  as  a  boundary  condition  which 
relates  the  mass  transfer  to  the  velocity  of  the  moving  surface  r  as  follows, 


m  =  p  r  (16). 

w  w 

In  the  case  of  an  isotropic  solid  in  which  the  thermal  properties  are  not  spatially 
dependent  but  arc  functions  of  temperature  only,  the  heat  flux  vector  Q  may  be  written 
in  terms  of  the  Fourier  relation 

Q  =  -  KVT  (17). 

This  relationship  then  gives  the  familiar  equation  of  heat  conduction 

pcp|-  =  V-  (KVr)  (18). 

In  this  equation  the  thermal  conductivity  K  is  a  scalar  temperature-dependent 
function. 


However,  when  one  considers  a  solid  which  is  anisotropic,  then  the  thermal  con¬ 
ductivity  is  position  dependent  as  well  as  temperature  dependent,  and  consequently  is 
a  second  order  tensor  having  nine  components.  In  the  particular  case  of  an  orthotro¬ 
pic  solid,  three  scalar  thermal  conductivity  coefficients  appear  and  the  Fourier  ex¬ 
pression  must  be  modified  as  follows  (4): 


-Q 


K. 


OT  -  JT_  +- 

9U!  "2  h2  du2  "3  h3 


(19) 


The  quantities  Kj,  K2,  K.}  are  the  three  mutually  perpendicular  thermal  conducti¬ 
vities;  hi,  hn,  I13  are  the  scale  factors  in  the  expression  for  the  increment  of  arc 
length  ds. 

9  09  2  2  2  2 

ds  =h1“du1  +h2  du2  +h3  du3  (20). 

For  specific  coordinate  systems,  the  foregoing  equations  may  be  reduced  to  the  results 
given  by  Nolan  and  Scala  (4).  For  example,  the  spherical  coordinate  system  appropriate 
to  a  nose  cap  region  may  be  written  for  the  case  of  axially-symmetric  flow: 
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CST  _  1  a  r„  2  st!  .  1  3  L  .  ,  3t1 

pcpl  1  IT  LKrr  V  +  1~~7  —  Ktsm*lJ  (21) 

r  J  r  sin  ip  T  L  T  T  -J 

At  the  stagnation  point  where  i|)  =  0,  it  is  to  be  noted  that  an  indeterminate  form  appeal's 
in  the  second  term  on  the  right  hand  side  of  Eq.  (21).  However,  upon  applying  L'Hopital's 
rule  to  this  term  one  finds  that  in  the  limit  as  ip  approaches  zero,  the  tern  is  finite  and 
has  the  value 

_2 _ rK  JT_~ 

2  dip  ip  dip  _ 

Hence,  at  the  stagnation  point  where  ip  =  0,  Eq.  (21)  may  be  written: 

5T  1  a  2  3T~I  ,  2  a  r  9T  ~| 

PCP  3t  r2  9r  15  r  SrJ+  2  dip  d<l>  J  22 

In  cylindrical  coordinates  Eq.  (15)  becomes: 

ar  i  a  L  arl  .  1  a  L  st  1  .  a  L  arl 


K  r  dT~ 

+  1  3  r, 

dT  1 

d 

K 

r  dr 

2  dip 

r  ^ 

Ip  d$_ 

dz 

z  3>z  _ 

If  one  assumes  that  the  external  flow  over  a  wing  or  fin  is  two-dimensional  so  that  cross 
flow  does  not  occur  along  the  span,  then  5T  may  be  set  equal  to  zero  in  Eq.  (23). 

dz 

We  may  now  discuss  boundary  conditions  appropriate  to  Eqs.  (21)  and  (23).  At 
time  t  =  0,  the  temperature  distribution  within  the  solid  may  be  either  constant  or  a 
prescribed  function  of  space.  That  is, 

T(r,  ip,  0)  =  f(r,  iP)  (24). 

But  because  of  chemical  reaction,  i.e. ,  oxidation  at  the  surface,  the  location  of 
the  outer  surface  is  time-dependent,  since  the  mass  loss  mw  is  a  function  of  surface 
position  for  given  temperatures  and  environmental  conditions.  In  the  case  of  the  hemi¬ 
sphere,  the  instantaneous  surface  recession  rate  may  then  be  expressed  by 


r(r,  ip,  t) 


(r,  ip,  t) 


Since  the  heat  flux  is  prescribed  at  the  surface,  Neumann  type  boundary  conditions 
are  to  be  imposed.  For  the  hemisphere  one  obtains: 


r  h  *■ l)]  “  "k~  K. 


4  4 

+  ae  T  -aeT 
t)  g  g  w  w. 


In  this  equation,  the  symbol  s(r,  ip,  t)  denotes  the  instantaneous  surface  position  of  the 
outer  boundary. 
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As  has  been  shown  by  Gibson  (31),  the  closed  form  solutions  for  moving  boundaries, 
even  in  one  dimension,  (Stefan's  problem),  are  restricted  to  a  very  few  cases  with 
Neumann  boundary  conditions  involving  a  time  dependent  heat  flux.  Consequently,  the 
foregoing  equations  with  appropriate  flight  parameters  have  been  programmed  for  high 
speed  digital  computer  solution.  Recently,  the  foregoing  transient  equations  were  ex¬ 
tended  at  General  Electric  to  include  the  case  of  multi-dimensional  heat  conduction  in 
conical  bodies  with  surface  oxidation  effects. 

THERMOSTRUCTURAL  RESPONSE 

Solutions  to  general  three-dimensional  problems  of  thermoelasticity  must  satisfy 
the  familiar  equations  of  equilibrium  and  conditions  of  compatibility.  Additionally, 
they  must  satisfy  the  prescribed  boundary  conditions  at  the  surface  of  the  body  under 
consideration.  These  equations  in  tensor  form  are  given  in  Ref.  (32)  and  are  listed 
here: 


a..  .  +  F.  =  0 
ij,  J  i 

€ij,  k-t  ek-t,  ij  fik,  j l  €j l,  ik  =  0  (27) 

o.A  .  =  T. 

U  i  i 

where  F.,  are  the  body  force  components,  T.  are  the  surface  force  components,  and 
are  the  normal  direction  cosines  of  the  outer  surface  in  the  x.  coordinate  directions. 

l 

These  equations  arc  based  on  purely  mechanical  and  geometrical  considerations 
and  therefore  would  apply  generally  to  any  elastic  solid,  whether  isotropic  or  aniso¬ 
tropic.  The  difference  that  exists  between  anisotropic  and  isotropic  materials  occurs 
in  the  stress-strain  relations.  For  the  most  general  case  of  an  anisotropic  elastic 
solid,  the  generalized  Hooke's  law  relation  contains  36  elastic  constants.  However,  it 
can  be  shown  by  arguments  based  on  thermodynamics  that  not  all  of  these  are  independent 
and  that  there  arc  21  independent  elastic  constants  for  the  general  case  of  anisotropy. 
When  thermal  effects  are  taken  into  ace  omit,  the  Hooke's  law  relations  contain  six  addi¬ 
tional  independent  coefficients  of  thermal  expansion.  If  the  solid  body  under  considera¬ 
tion  exhibits  structural  symmetry  with  respect  to  a  plane  or  about  an  axis,  the  number 
of  independent  elastic  and  thermal  constants  is  reduced  considerably. 

Pyrolvtieally  deposited  materials  exhibit  symmetry  with  respect  to  three  orthogonal 
planes  and  the  c  direction  axis  which  is  normal  to  the  ab  plane  of  deposition.  Love  (33) 
refers  to  a  material  exhibiting  such  symmetry  as  "transversely  isotropic".  However, 
others  have  suggested  the  terms  "transversely  anisotropic"  and  "monotropic"  as  more 
descriptive.  It  is  shown  in  Ref.  (32)  that  a  material  which  exhibits  symmetry  with 
respect  to  three  orthogonal  planes  is  called  an  orthotropic  material  and  has  the  follow¬ 
ing  stress-strain  relations: 
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Talcing  the  ab  plane  as  the  plane  of  deposition,  Eq.  (28)  can  be  further  specialized  for 
the  monotropic  case  by  requiring  elastic  symmetry  with  respect  to  the  c  axis.  Thus, 
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(30). 


It  is  more  convenient  to  write  the  monotropic  stress-strain  relations  in  terms  of  pro¬ 
perties  in  the  plane  of  deposition;  E,  v,  and  a  (for  convenience  the  subscript  a  is  no 
longer  used),  and  anisotropy  ratios  jS,  /i,  k,  and  X  are  defined  by: 


=  _G__  =  E  1 

M  G  1  +  v  G 
c  c 

(31). 

K  =  — — — 

0! 

C 


V 


C 


Note  that  for  an  isotropic  material  the  anisotropy  ratios  become  equal  to  unity.  Thus, 
the  thermoelastic  stress-strain  equations  for  a  monotropic  material  with  the  c  axis 
as  the  axis  of  symmetry  can  be  expressed  in  the  following  form: 
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Many  of  the  geometries  that  are  of  practical  interest  in  aerospace  structural  design 
have  dimensions  in  the  c  direction  which  are  usually  small  compared  to  those  in  the 
ab  plane,  such  as  leading  edge  and  nose  can  configurations.  For  these,  the  assumption 
of  plane  stress  and  plane  strain  can  be  made  which  will  considerably  simplify  any 
analysis. 


Consider  a  long  hollow  cylindrical  shell  of  circular  cross  section  subjected  to  a 
temperature  distribution  which  is  a  function  of  radius  only,  see  Fig.  (18).  The  cylin¬ 
drical  coordinates  r,  0,  z,  coincide,  respectively,  with  the  c,  a,  b  coordinates  of  the 
monotropic  material.  Because  of  the  radial  symmetry  of  the  geometry  and  temperature 
distribution,  the  assumption  of  plane  strain  for  interior  cross  sections  may  be  made. 
The  corresponding  equations  of  equilibrium,  compatibility,  and  stress  strain  are: 
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The  above  equations  can  be  reduced  to  one  governing  equation  for  the  radial  stress: 
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Fig.  18.  Cylindrical  Shell  Coordinate  Syste 


The  tangential  stresses  can  be  computed  directly  from  the  equation  of  equilibrium,  Eq. 
(33).  In  certain  cases  it  may  be  possible  to  ignore  the  valuation  of  mechanical  proper¬ 
ties  with  temperature.  For  such  a  case,  Eq.  (34)  reduces  to  a  much  simpler  form: 
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For  the  isotropic  cylinder,  the  anisotropy  ratios  )3,  X,  and  K  arc  equal  to  one,  and  Eq. 
(34)  reduces  to: 
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For  constant  mechanical  and  thermal  properties,  Eq.  (37)  becomes: 
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The  most  significant  difference  between  the  behavior  of  the  monotropic  cylinder,  Eqs. 
(34)  mid  (35),  and  the  isotropic  cylinder,  Eqs.  (37)  mid  (38),  is  associated  with  the 
fact  that  the  monotropic  equations  possess  a  term  on  the  right-hand  side  of  the  differ¬ 
ential  equation  which  is  proportional  to  the  local  temperature,  in  addition  to  a  term 
involving  a  variation  of  temperature  with  radius.  A  monotropic  cylinder  free  from  any 
external  constraint  will  then  develop  thermal  stresses  while  moving  from  one  uniform 
temperature  environment  to  another.  This  is  caused  by  the  unequal  thermal  expansion 
coefficients  a  mid  a  which  give  rise  to  internal  geometric  constraints  causing  thermal 
stresses  to  be  created.  Another  interesting  feature  of  the  monotropic  body  is  that  the 
nature  of  the  stress  state  depends  on  the  anisotropy  ratios  /?,  X,  and  k.  Thus,  when 
a  hollow  cylinder  is  subjected  to  a  temperature  gradient  through  the  thickness  such  that 
the  outer  radius  is  at  a  higher  temperature  than  the  inner  radius,  it  is  usually  antici¬ 
pated  that  compressive  tangential  stresses  will  exist  at  the  outer  radius  and  tensile 
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tangential  stresses  at  the  inner  radius.  The  corresponding  radial  stress  distribution 
would  then  be  tension  throughout  the  thickness.  This  stress  distribution  always  occurs 
for  the  isotropic  cylinder.  However,  for  the  monotropic  cylinder,  depending  on  the 
geometry  and  anisotropy  ratios,  a  complete  reversal  can  occur  —  that  is,  the  radial 
stress  distribution  may  be  compressive  throughout  the  thickness  and  the  tangential 
stresses  at  the  inner  and  outer  surfaces  are  compressive  and  tensile,  respectively. 
This  is  an  extremely  important  factor  in  design,  since  the  compressive  c  direction 
strength  of  pyrolytic  graphite  materials  is  much  greater  than  the  tensile  c  direction 
values. 

THERMAL  STRESSES  IN  A  SPHERICAL  SHELL 


Consider  now  a  spherical  shell  of  inner  radius  Rj  and  R  subjected  to  a  tempera¬ 
ture  distribution  which  is  a  function  of  radius  only,  Fig.  (19).  The  deposition  geometry 
is  such  that  the  middle  plane  of  the  shell  coincides  with  the  plane  of  deposition.  Thus, 
the  governing  equations  of  equilibrium,  compatibility,  and  stress  strain  are: 
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The  above  equations  can  be  combined  in  the  same  manner  as  that  used  for  the  cylinder, 
resulting  in  the  following  differential  equation  for  the  radial  stresses  in  a  monotropic 
sphere: 
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The  tangential  stresses  can  be  computed  directly  from  Eq.  (39).  If  the  mechanical  pro¬ 
perties  are  held  constant  with  temperature,  Eq.  (40)  reduces  to: 
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The  isotropic  forms  of  Eqs.  (40)  and  (41)  are 
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Inspection  of  Eqs.  (40)  to  (45)  indicates  that  the  comments  made  earlier  regarding  the 
differences  in  behavior  between  the  monotropic  and  isotropic  cylinder  apply  equally  well 
to  the  spherical  shell. 


RESIDUAL  STRESS 

Among  die  many  unusual  characteristics  which  arise  from  the  monotropic  nature 
of  pyrolytically  deposited  materials  is  the  extreme  importance  of  residual  stresses  in 
producing  sound  "as  deposited"  hardware.  It  is  apparent  that  the  thermal  anisotropy, 
which  is  desirable  from  a  thermal  protection  viewpoint,  can  result  in  residual  stresses 
which  can  be  detrimental  to  the  structural  integrity  of  the  component.  That  is,  in  addi¬ 
tion  to  the  primary  concern  of  affecting  structural  integrity  during  extreme  re-entry 
environmental  loadings,  these  residual  stresses  present  a  fundamental  problem  insofar 
as  their  effect  upon  the  production  yield  of  flight  hardware,  including  producing  "as 
deposited"  shapes,  machining  to  flight  thickness,  and  machining  of  holes.  The  effects 
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which  contribute  to  the  residual  stress  state  of  structural  elements  fabricated  from  py¬ 
rolytic  materials  are  listed  below. 


1.  General  Anisotropy 

2.  Growth 

3.  Thickness  Temperature  Drop 

4.  Nodules 

5.  Radial  Anisotropy 

6.  Non-Uniform  Axial  Temperature 

7.  Asymmetric  Effects 

8.  Geometric  Discontinuities 

General  anisotropy  is  associated  with  those  residual  stresses  which  arise  upon 
cooling  to  room  temperature  from  deposition  temperature  (about  4000°F).  This  effect 
is  one  of  the  major  causes  of  residual  stresses  in  pyrolytic  materials  and  is  due  to  the 
monotropic  character  of  the  material.  On  thin-walled  pyrolytic  graphite  nose  caps,  the 
circumferential  and  meridional  stresses  arc  compressive  on  the  outside  and  tensile 
on  the  inside.  Radial  stresses  are  tensile  and  reach  their  maximum  value  in  the  in¬ 
terior  of  the  thickness.  On  cylindrical  and  conical  flare  sections,  the  circumferential 
stresses  are  compressive  on  the  outside  surface  and  tensile  on  the  inside.  Meridional 
stresses  are  smaller  in  magnitude  and  of  opposite  sign.  The  radial  stresses  are  of 
the  same  sense  as  those  in  the  nose  cap  region.  Since  the  allowable  tensile  strengths 
of  pyrolytic  graphite  are  low  in  the  c  direction,  any  effect  which  tends  to  delaminate 
the  material,  such  as  general  anisotropy,  must  be  carefully  studied. 

The  biaxial  stresses  which  result  from  the  anisotropic  expansion  coefficient  and 
the  high  deposition  temperature  are  theoretically  calculable  for  geometrically  homo¬ 
geneous  shapes,  such  as  cylinders  and  spheres,  from  the  equations  presented  earlier. 
For  a  given  geometry,  the  general  anisotropy  stresses  are  a  function  of  the  thickness  to 
radius  ratio.  For  a  sphere,  (lie  maximum  tensile  and  compressive  tangential 
stresses  increase  linearly  with  t/R,  Fig.  (20).  The  maximum  radial  stress  in  the 
sphere  increases  approximately  as  the  square  of  t/R,  Fig.  (21). 

The  major  contribution  to  residual  stress,  other  than  pure  anisotropy,  is  the  stress 
induced  by  variable  lattice  transformation  of  "growth"  during  the  deposition  process. 

In  order  to  explain  these  stresses,  it  is  first  necessary  to  discuss  lattice  transforma¬ 
tion  and  the  characteristics  of  the  deposition  process.  Lattice  transformation  is  a  phe¬ 
nomenon  of  dimensional  change  of  the  lattice  due  to  soaking  at  high  temperature.  The 
dimensional  change  in  the  macroscopic  piece  is  therefore  an  increase  in  length  and 
circumference  and  a  decrease  in  thickness.  For  any  given  temperature,  there  is  n 
definite  maximum  change  or  transformation.  As  the  temperature  increases,  the  amount 
of  change  increases.  The  kinetics  of  transformation  are  first  order  with  respect  to  tem¬ 
perature,  so  that,  as  the  temperature  increases,  the  rate  of  transformation  increases. 

A  plot  of  dimensional  change  versus  time  will  show  an  exponential  rise  to  an  asymptotic 
value  for  any  given  temperature. 
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The  time  of  deposition  for  typical  pyrolytic  graphite  hardware  can  vary  from  as  low 
as  10  to  12  hr.  to  as  high  as  130  hr.  Considering  any  given  run,  the  initial  layer  is 
deposited  within  several  micro-seconds;  the  succeeding  layers  are  built  up  continuously 
at  an  essentially  constant  rate  until  the  last  layer  is  deposited  and  then  the  furnace  is 
turned  off.  Immediate  cooling  takes  place.  The  first  layers  are  exposed  to  the  deposi¬ 
tion  temperature  for  the  total  time,  and  each  succeeding  layer  is  exposed  for  less  time 
until  the  final  layer  which  is  essentially  unexposed.  The  time  temperature  history  of 
each  layer  is,  therefore,  different  and  the  different  "growths"  set  up  biaxial  distribu¬ 
tions  which  are  opposite  in  sign  to  the  general  anisotropy  residual  stresses,  which  is 
an  important  consideration  in  designing  sound  "as  deposited"  hardware. 

A  further  factor  to  consider  here  is  the  role  of  the  mandrel  in  the  "growth"  stress 
picture.  The  mandrel,  being  commercial  graphite,  does  not  undergo  dimensional  change 
at  the  deposition  temperature  and  therefore  acts  as  a  constraining  fixture.  The  growth 
causes  the  mandrel  and  any  previously  deposited  material  to  be  expanded.  When  the 
deposition  layer  is  sufficiently  thick,  sufficient  tension  is  built  up  in  the  manorel  to 
fracture  it.  The  mandrel  falls  off  and  the  subsequent  deposition  takes  place  on  the  pre¬ 
viously  deposited  material.  As  deposition  continues,  the  first  deposited  layers  are 
pushed  farther  and  farther  into  tension  by  the  outer  layers.  The  final  stress  distribu¬ 
tion,  as  a  result  of  growth,  is  a  residual  tensile  stress  on  the  mandrel  side.  The  mag¬ 
nitude  of  the  compressive  stress  at  the  last  deposited  side  depends  on  the  deposition 
rate  in  addition  to  being  a  characteristic  of  the  material,  and  is  of  the  order  of 
10,000  psi. 

Because  of  (lie  complexity  of  the  growth  mechanisms,  these  residual  stresses  can¬ 
not  be  computed  as  simply  as  those  arising  from  general  anisotropy.  The  most  conve¬ 
nient  procedure  is  to  obtain  a  growth  stress  formulation  empirically.  This  procedure 
is  dependent,  however,  upon  the  complete  cognizance  of  all  factors  affecting  growth 
stresses.  Levy  (34)  lias  proposed  a  growth  stress  theory  which  assumes  a  very  rapid 
rate  of  growth  that  sets  up  equivalent  pressures  and  hence  stresses  as  a  result  of  the 
inherent  restraints  agains  this  growth.  In  order  to  obtain  a  numerical  solution  using 
this  formulation,  it  is  necessary  to  know  the  ultimate  strength  of  the  mandrel  material 
at  the  temperature  at  which  it  breaks,  the  mandrel  thickness,  and  the  stress  in  the 
plane  of  a  layer  due  to  growth  without  changes  in  dimensions,  i.e.,  fully  constrained 
growth  stress  —  in  addition  to  the  geometry  of  the  component.  Further  work  in  this 
area  is  currently  understudy  at  the  General  Electric  Re-Entry  Systems  Department. 

Since  a  vapor-deposited  material  such  as  pyrolytic  graphite  has  good  insulation 
characteristics  in  the  e  direction,  each  successive  layer  is  deposited  at  a  slightly  lower 
temperature  when  a  constant  externally  applied  heat  source  is  used.  Accordingly,  after 
the  entire  thickness  has  been  deposited,  for  computational  purposes,  the  shell  can  be 
construed  as  having  a  temperature  drop  across  the  thickness.  The  magnitude  of  this 
temperature  drop  would  be  proportional  to  the  thickness  deposited.  However,  the 
magnitude  of  this  drop  for  the  thicknesses  of  structures  under  consideration  is  not 
large,  and  its  effect  upon  the  residual  stress  picture  is  not  a  major  one.  Analytical 
techniques  are  available  and  can  be  easily  used  to  compute  these  effects,  which  arise 
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as  the  shell  cools  down  from  this  temperature  distribution  to  an  isothermal  room  tem¬ 
perature  condition. 

Nodules  (abnormal  growths  in  pyrolytic  graphite  that  disrupt  the  continuity  of  the 
parallel  deposit  planes)  were  a  problem  during  the  early  stages  of  development.  They 
are  a  localized  source  of  residual  stress,  as  well  as  a  point  of  weakness,  in  many 
cases  caused  by  the  presence  of  a  soot  particle  in  the  How.  The  analytical  treatment 
by  Levy  (34),  verified  by  experiment,  showed  that  these  stresses  —  although  serious 
for  comparatively  large  nodules  --  become  secondary  as  the  inclusion  size  diminishes. 
As  the  causes  of  their  formation  became  known,  methods  were  introduced  in  the  depo¬ 
sition  process  to  minimize  (and,  in  some  cases,  virtually  eliminate)  them  from  the 
finished  product.  This  was  accomplished  by  the  use  of  a  better-grade  graphite  for  man¬ 
drels  and  other  parts,  control  of  mandrel  surface  roughness,  and  thorough  removal  of 
surface  dust  from  all  parts. 

It  has  also  been  determined  experimentally  that  the  material  properties  sometimes 
vary  through  the  thickness  of  the  deposit.  Thus,  strictly  speaking,  the  material  is 
transversely  nonhomogeneous.  The  stresses  that  result  from  this  source  upon  cooling 
to  room  temperature  are  comparatively  small,  the  principal  effect  being  a  small  change 
in  the  radial  stress.  In  addition,  since  the  gas  injector  is  a  heat  sink,  and  since  the 
radiative  heat  loss  will  vary  inversely  as  the  distance  to  the  heat  sink,  a  nonuniform 
axial  deposition  temperature  exists.  Again,  the  magnitude  of  this  temperature  varia¬ 
tion  is  not  large  and  has  a  small  effect  upon  tlie  axial  residual  stresses.  Another 
source  of  residual  stress  is  associated  with  asymmetric  effects  such  as  circumferential 
variations  of  thickness,  temperature  of  deposition,  and  material  properties.  Because 
of  the  extremely  random  nature  of  these  variations,  the  residual  stresses  arising  from 
sources  cannot  be  determined  analytically  with  any  accuracy.  However,  by  proper 
process  control  during  deposition  these  variations  can  be  minimized  and  have  been 
reduced  in  importance. 

Deposition  of  pyrolytic  materials  on  mandrels  in  the  form  of  shells  of  revolution 
with  finite  lengths  can  complicate  the  residual  stress  pattern  due  to  discontinuities  asso¬ 
ciated  with  sudden  changes  in  shell  thickness,  juncture  of  two  shells  of  revolution,  and 
a  free  edge  at  the  end  of  a  shell.  While  the  structural  effects  associated  with  these 
discontinuities  as  well  as  those  caused  by  thermal  gradients  tend  to  be  absorbed  via 
an  essentially  plane  strain-stress  action  in  the  plane  of  deposition  one  cannot  utilize 
the  conventional  theory  of  shells.  Additional  important  effects  caused  by  the  aniso¬ 
tropy  of  the  thermal  expansion  coefficients,  and  the  low  shear  rigidity  between  layers 
must  be  taken  into  account.  Use  of  the  classical  isotropic  shell  equations  can  result 
in  serious  errors  for  may  cases  of  practical  interest.  The  work  at  GE  MSD  asso¬ 
ciated  with  the  development  of  a  monotropic  shell  theory  will  be  discussed  later. 


General  anisotropy,  growth,  and  geometric  discontinuities  are  the  major  causes 
of  residual  stresses  in  "as  depositied"  hardware.  These  effects  were  not  well  under¬ 
stood  during  the  early  stages  of  pyrolytic  material  development,  and  as  cracks  and 
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delaminations  began  to  occur  in  production  material  a  considerable  effort  was  spent  on 
the  study  of  the  problem.  Consequently,  manufacturing  limits  associated  with  minimum 
curvature,  size,  and  thickness  were  established  and  a  series  of  nondestructive  tests 
to  evaluate  the  soundness  of  the  material  were  developed.  These  nondestructive  tests 
include  dyecheck,  videogauging,  imagescope,  and  other  ultrasonic  and  x-ray  methods. 
The  state  of  the  art  has  developed  to  the  point  where  today's  manufacturers  can  deliver 
materials  of  thicknesses  considered  impossible  several  years  ago. 

In  order  to  further  increase  the  limiting  thickness  radius  ratios  for  pyrolytic  ma¬ 
terial  components,  studies  have  been  conducted  at  GE-RSD  to  investigate  the  interaction 
between  growth  and  general  anisotropy  stresses.  This  would  apply  to  those  portions 
which  are  beyond  die  effects  of  free  edges  or  similar  discontinuities,  such  as  the  stag¬ 
nation  region  of  nose  caps  and  leading  edges.  Since  growth  stresses  are  opposite  in 
sign  to  die  "cool  down"  stresses,  calculations  were  made  to  define  die  mandrel  tech¬ 
nology  necessary  to  develop  compressive  radial  growtii  stresses  of  sufficient  magnitude 
to  overcome  the  tensile  radial  general  anisotropy  stresses  and  result  in  a  net  radial 
compressive  stress  state.  This  would  ensure  a  residual  stress  state  that  would  be 
free  of  any  tendency  for  the  material  to  delaminate.  The  calculations  indicate  that  it 
is  indeed  possible,  for  certain  configurations  of  practical  interest,  to  specify  a  man¬ 
drel  geometry  diat  will  ensure  net  compressive  radial  stresses,  tangential  stresses 
within  the  "a"  direction  allowables  of  the  material,  and  an  overall  residual  stress  state 
which  is  opposite  to  the  environmental  stress  state.  Work  is  currently  under  way 
to  verify  experimentally  die  dieorctical  calculations.  This  interaction  between  growdi 
and  general  anisotropy  is  shown  in  Figs.  (22)  and  (23)  for  a  6  in.  boron  pyrolytic  graphite 
sphere  of  0.  3  in.  diick.  The  discontinuity  in  the  curves  is  due  to  the  fracturing  of  die 
mandrel  which  occurs  when  the  deposit  is  0.05  in.  thick.  The  growth  stresses  result 
in  a  50  per  cent  reduction  in  die  peak  radial  tensile  stress. 

DESIGN  ANALYSIS  TECHNIQUES 

In  order  to  properly  design  practical  structural  components  fabricated  out  of  pyro¬ 
lytic  materials,  it  is  necessary  to  have  available  design  analysis  techniques  which 
reliably  predict  the  residual  and  environmental  stress  conditions  in  the  component. 
Limiting  the  discussion  to  shell  structures,  the  governing  factors  are  dien  the  spatial 
variation  of  the  temperature  distribution  and  the  geometry  of  the  component.  Assum¬ 
ing  that  diesc  can  be  expressed  in  some  mathematical  form,  then  it  is  indeed  possible 
to  cast  the  thermal  stress  problem  for  monotropic  structures  in  a  mathematical  form 
that  represents  an  exact  three-dimensional  approach  to  the  problem.  The  major  pro¬ 
blem  is  in  die  solution  of  the  differential  equations  which  will  yield  die  stresses  and 
displacements  of  the  structure.  Even  with  the  use  of  modern  computer  facilities, 
stress  analysts  have  been  unable  to  obtain  versatile  and  useful  design  analysis  tech¬ 
niques  for  isotropic  shell  structures  by  a  frontal  attack  on  the  exact  three-dimensional 
equations  of  isotropic  elasticity.  Consequently,  it  has  been  necessary  to  develop  a 
dieory  for  shells  in  which  certain  simplifications  have  been  made  in  order  to  make  the 
mathematical  difficulties  less  severe.  The  validity  of  these  simplifications  for  thin 
shells  have  been  so  sufficiently  verified  by  lest  tiiat  stress  analysts  today  use  the 
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Figo  23„  Net  Tangential  Residual  Stresses  in  a  Spherical  Shell 


engineering  theory  for  isotropic  shell  structures  without  hesitation. 

The  exact  three-dimensional  problem  for  monotropic  shell  structures  presents 
as  much,  if  not  more,  of  a  problem  than  the  isotropic  problem.  Only  in  certain  spe¬ 
cial  limited  cases  as  discussed  previously  can  one  obtain  exact  solutions.  Thus,  it  is 
necessary  to  develop  an  engineering  theory  for  monotropic  shells  that  will  enable  the 
stress  analyst  to  compute  the  thermal  and  environmental  stresses  in  practical  structures. 

A  theory  for  monotropic  shells  is  under  development  at  GE-RSD  and  the  initial  re¬ 
sults  are  very  encouraging.  McDonough  (35)  has  developed  a  shell  theory  for  monotropic 
shells  of  revolution  which  takes  into  account  the  unusual  thermostructural  characteris¬ 
tics  of  monotropic  materials.  The  theory  was  initially  developed  for  a  spherical  shell 
and  then  extended  to  a  general  shell  of  revolution.  The  complete  three-dimensional 
equations  of  elasticity  were  modified  using  appropriate  shell  approximations  and  the 
governing  shell  equations  were  obtained  by  direct  integration.  A  comparison  was  made 
between  the  predictions  of  the  shell  theory  and  the  complete  three-dimensional  solution 
for  a  symmetrically  heated  sphere,  Figs.  (24)  and  (25).  The  agreement  is  very  good. 

The  shell  technique  is  currently  being  programmed  for  digital  computation,  for  the 
case  of  the  general  shell  of  revolution. 

FUTURE  RESEARCH  AND  DEVELOPMENT 

The  attractiveness  of  pyrolytic  graphite  from  a  thermal  and  structural  viewpoint 
for  certain  hypersonic  vehicles  has  been  established,  see  Refs.  (3)  to  (G).  One  must 
not  imply,  however,  that  the  ultimate  has  been  reached  for  this  type  of  material;  for 
example,  evidence  has  already  been  obtained  (G),  (23),  that  alloying  the  pyrolytic  gra¬ 
phite  with  a  small  amount  of  boron  increases  both  the  reaction  rate  controlled  resis¬ 
tance  to  oxidation  and  the  tensile  strength  at  high  temperature.  The  possibility  of  pro¬ 
ducing  other  combinations  of  pyrolytic  graphite  and  refractory  substances  such  as 
hafnium,  tantalum,  etc.  has  already  been  investigated  (29)  at  the  General  Electric 
Company,  Missile  and  Space  Division,  as  well  as  by  others.  Each  combination  of 
elements  yields  certain  advantages  not  present  in  the  others.  For  each  of  these  newer 
nuUnrfalff,  il  uilll  In*  to  Uhr  miLtikfttm'cal  ynrL  rmance  in 

a  hypersonic  environment. 

The  advantages  of  varying  the  degree  of  anisotropy  in  pyrolytic  graphite  to  minimize 
radial  stresses  during  hypersonic  Right  have  been  discussed  before  (5)  in  a  theoretical 
way.  One  would  anticipate  that  the  ability  to  produce  this  type  of  hardware  would  be 
well  worth  the  effort  invested  and  would  open  up  possibilities  for  missions  not  hereto¬ 
fore  considered. 

An  additional  area  of  further  development  lies  in  the  technique  of  high  temperature 
bonding  pyrolytic  graphite  to  other  structural  members  or  itself.  Work  at  GE  MSD  has 
produced  useful  bonding  materials  and  techniques  which  have  been  successfully  tested 
at  temperatures  ranging  up  to  2500°F.  Since  pyrolytic  graphite  has  a  thermal  capa¬ 
bility  in  excess  of  this  temperature,  one  may  conclude  that  effort  in  the  art  of  bonding 
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has  progressed  substantially  over  that  which  existed  only  two  years  ago.  Present  indi¬ 
cations,  based  on  very  recent  effort  at  GE  MSD  (23),  are  that  bonds  for  much  higher 
temperatures  can  be  successfully  produced  and  applied  to  pyrolytic  graphite.  One 
should  also  comment  in  passing,  that  it  may  well  be  the  case  in  actual  practice  that  the 
bond  line  temperature  might  not  exceed  say  1000°F.  For  many  hypersonic  missions, 
sufficient  time  may  not  elapse  to  raise  the  bond  line  temperature  any  significant  amount 
at  all.  This  is  a  consequence  of  the  low  "c"  direction  thermal  conductivity  and  high 
oxidation  resistance  of  pyrolytic  materials,  (4). 

Future  effort  in  the  area  of  material  property  investigations  should  be  directed 
toward  categorizing  the  performance  data  with  a  sound  and  identifiable  material  struc¬ 
ture  (23).  For  instance,  it  has  already  been  established  that  the  coefficient  of  thermal 
expansion  is  substantially  affected  by  changes  in  the  material  structure  (23).  This 
property  plays  a  significant  role  in  determining  the  stress  level  both  during  manufac¬ 
ture  and  hypersonic  flight,  (6).  It  is  apparent  that  more  attention  must  be  paid  to  this 
problem  to  assure  one's  self  that  flight  hardware  is  indeed  representative  of  the  spe¬ 
cimens  used  to  measure  the  material  properties  upon  which  the  design  is  based. 

CONCLUSIONS 

During  the  past  several  years,  effort  in  the  development  of  pyrolytic  graphite  has 
produced  a  material  entirely  different  than  that  produced  in  the  early  days  of  1959-1960. 

Knowledge  regarding  the  manufacturing  techniques,  deposition  parameters,  man¬ 
drel  design  and  characterization  of  the  material  structure  has  been  significantly  in¬ 
creased  since  1961.  The  art  has  progressed  to  the  point  where  intricate  shapes  may 
now  be  made  in  varying  thicknesses,  whereas  just  one  year  ago  such  a  product  was 
thought  to  be  impossible  using  the  vapor  deposition  process. 

Reliable  mechanical  data  have  now  been  obtained  at  temperatures  to  5000°F.  These 
data  show  that  the  strength  at  high  temperatures  is  indeed  beyond  any  other  known 
material.  Analytical  techniques  not  known  two  years  ago  have  now  been  developed  and 
confirmed  to  predict  not  only  the  hypersonic  flight  response  but  the  stress  state  of  the 
manufactured  piece  as  well.  As  a  consequence  of  these  techniques,  an  awareness  of 
the  importance  of  relating  a  well-characterized  structure  to  the  material  properties 
has  resulted,  and  renewed  effort  is  already  underway  to  correlate  this  data. 

Even  though  a  complete  array  of  high  temperature  data  is  not  available,  enough 
data  now  exists  to  enable  one  to  design  experimental  flight  hardware  using  pyrolytic 
graphite. 


SYMBOLS 

c  specific  heat  at  constant  pressure 
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activation  energy,  modulus  of  elasticity 

body  force 

shear  modulus 

enthalpy 

scale  lengths 

stagnation  enthalpy 
specific  reactivity 
effective  collision  frequency 

thermal  conductivity 
direction  cosines 

mass  transfer  rate 
order  of  the  reaction 
local  static  pressure 
heat  transfer  rate 
radius,  radial  component 
recession  rate 
gas  constant 
nose  radius 

instantaneous  location  of  the  moving  boundary 
time 

temperature,  surface  force 

unit  vector 

coordinate 

mole  fraction 

velocity 

axial  component 
defined  in  Eq.  (31) 
emissivity,  strain 
angular  component 
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K 

defined  in  Eq.  (31) 

X 

defined  in  Eq.  (31) 

A 

angle  of  yaw 

M 

defined  in  Eq.  (31) 

V 

Poisson's  ratio 

P 

density 

a 

stress 

$ 

defined  in  Eq.  (36) 

1 1) 

meridional  angle 

'I' 

defined  in  Eq.  (36) 

Cl 

defined  in  Eq.  (34) 

Subscripts 

a 

a-b  plane 

D 

diffusion  controlled  regime 

e 

edge  of  boundary  layer 

B 

gas  cap 

i 

species  i 

ij,  i a 

tensor  components 

R 

reaction  rate  controlled  regime 

S 

sublimation  regime 

\V 

wall 
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summary^.  -s  pUrpose  0f  fhfg  paper  to  review  the  recent 
history  of  the  development  of  pyrolytic  graphite  and  its 
thermal3  aerothermochemical,  and  structural  properties 
with  regard  to  its  application  to  hypersonic  vehicles. 

Pyrolytic  graphite  is  not  really  new  since  its  dis¬ 
covery  dates  back  to  Thomas  Edison,  however,  this  class 
K)f  refractory  materials  has  been  "rediscovered"  recently  - 
because  it  possesses  certain  desirable  properties  at  ex¬ 
tremely  high  operating  temperatures,  including  a  high 
strength  to  weight  ratio,  low  oxidation  rate,  and  a  low 
transverse  coefficient  of  thermal  conductivity. 

Consequently,  a  thorough  investigation  of  the  prop¬ 
erties  of  pyrolytic  graphite  is  being  carried  out.  In  ad¬ 
dition,  theoretical  models  are  being  developed  so  that  its 
aerothermochemical  and  aer othermostructural  behavior 
I  can  be  predicted  with  confidence  for  arbitrary  hypersonic 
missions.  In  this  paper,  the  currently  available  property 
data,  and  analytical  and  numerical  techniques  are  present¬ 
ed  and  discussed  critically,  and  areas  for  future  research 
and  development  are  delineated. 
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